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Place cells in the rat hippocampus play a key role in creating the animal’s internal representation 
of the world. During active navigation, these cells spike only in discrete locations, together encoding 
a map of the environment. Electrophysiological recordings have shown that the animal can revisit 
this map mentally, during both sleep and awake states, reactivating the place cells that fired during its 
exploration in the same sequence they were originally activated. Although consistency of place cell 
activity during active navigation is arguably enforced by sensory and proprioceptive inputs, it remains 
unclear how a consistent representation of space can be maintained during spontaneous replay. We 
propose a model that can account for this phenomenon and suggests that a spatially consistent replay 
requires a number of constraints on the hippocampal network that affect its synaptic architecture and 
the statistics of synaptic connection strengths. 
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I. INTRODUCTION 

In the course of learning a spatial environment, an animal forms an internal representation of space 
that enables spatial navigation and planning (T). The hippocampus plays a key role in producing this map 
through the activity of location-specific place cells |@. At the neurophysiological level, these place cells 
exhibit spatially selective spiking activity. As the animal navigates its environment, the place cell fires only 
at a discrete location—its place field (Figure |T]A-B ). It is believed that the entire ensemble of place cells 
serves as a neuronal basis of the animal’s spatial awareness 0|4[- 

Remarkably, place cells spike not only during active navigation but also during quiescent wake states 
12 |6l and even during sleep 00. For example, the animal can “replay” place cells in sequences that 
correspond to the physical routes traversed during active navigation Ifl0l - fl2l or “preplay” sequences that 
represent possible future trajectories, either in direct or reversed order, while pausing at a decision point lfl4l 
mi. This phenomenon implies that, after learning, the animal can explore and retrieve spatial information 
by cuing the hippocampal network IIT61IT71 . which may in turn be viewed as a physiological correlate of 
“mental exploration” Ifl8l [T9l . 

It bears noting, however, that the actual functional units for spatial information processing in the hip¬ 
pocampal network are not individual cells but repeatedly activated groups of place cells known as cell 
assemblies (see |j20ll and Figure [Ip). Although the physiological properties of the place cell assemblies 
remain largely unknown, it is believed that the cells constituting an assembly synaptically drive a certain 
readout unit downstream from the hippocampus. In the “reader-centric” view, this readout neuron—a small 
network or, most likely, a single neuron—is what actually defines the cell assembly, by actualizing the in¬ 
formation provided by its activity Il20ll . The identity of the readout neurons in some cases is suggested by 
the network’s anatomy. For example, there are direct many-to-one projections from the CA3 region of the 
hippocampus to the CA1 region. Since replays are believed to be initiated in CA3 lfl3llT4l . this implies that 
the CA1 place cells may serve as the readout neurons for the activity of the CA3 place cells. Assuming that 
contemporaneous spiking of place cells implies overlap of their respective place fields (Figure [jp-B), it is 
possible to decode the rat’s current location from the ongoing spiking activity of a mere 40-50 neurons liTTll . 
This suggests that the readout neurons may be wired to encode spatial connectivity between place fields by 
responding to place cell coactivity (see Figure [1}\-C and I l22lf24l ). 

A natural assumption underlying both the trajectory reconstructing algorithms lf2Ti and various path 
integration models 0 [27] [29]] is that the representation of spatial locations during physical navigation is 
reproducible. If the rat begins locomotion at a certain location and at a certain moment of time, to, and then 
returns to the same location at a later time, t\, then the population activity of the place cells at to and t\ is the 
same. Similarly, if spatial information is consistently represented during replays, then the activity packet in 
the hippocampal network should be restored upon “replaying” a closed path. Whereas the correspondence 
between place cell activity and spatial locations (i.e., place fields) during physical navigation is enforced 
by sensory and proprioceptive inputs E71 . the consistency of spatial representation during replay must be 
attributable solely to the network’s internal dynamics J28J. 

Here we develop a model that accounts for how a neuronal network could maintain consistency of 
spatial information over the course of multiple replays or preplays. This model is based on the discrete 
differential geometry theory developed in ll30Tl . which reveals that key geometric concepts can be expressed 
in purely combinatoric terms. The choice of this theory is driven in part by recent work that indicates 
that the hippocampus provides a topological framework for spatial information rather than a geometric or 
Cartesian map Il24l426ll . 

The results suggest that to maintain consistency of spatial information during path replay, the synaptic 
connections between the place cells and the readout neurons must adhere to a zero holonomy principle. 
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FIG. 1 : Place cells, place fields and cell assemblies. A. A schematic representation of the spike trains produced by 
three place cells, C\, C2 and C3. The two red rectangles mark the periods during which the cells are coactive ESI- ». 
The gold, green and blue areas represent place fields. Place cell firing rate is maximal at the center of the place field 
and attenuates towards its periphery; this pattern can be closely approximated by Gaussian distribution. Place cell 
cofiring reflects overlap between respective place fields: cells c 1 and ci are coactive in the location Z12, cells C2 and 
C3 are coactive in the location I23 and so on. The red links mark distances between the centers of the place fields and 
the triple overlap domain, /123 (dark region in the center). C. A schematic representation of a cell assembly: the three 
place cells on the top synapse onto a readout neuron (red dot), which activates within the cell assembly field Z123. D. 
During replay, the place cells repeat on a millisecond time scale the order of spiking that they exhibit during active 
navigation. 


II. THE MODEL 


The simplicial model of the cell assembly network. A convenient framework for representing a pop¬ 
ulation of place cell assemblies is provided by simplicial topology 132I4341 . In this approach, an assembly 
of d + 1 place cells, {c; 0 , c,,,..., c/ rf ), is represented by a d-dimensional abstract simplex (not to be confused 
with a geometric simplex) containing d + 1 vertexes, cr = [v, 0 , v/,,..., v/J, where each vertex, v;, corresponds 
to a place cell c,- (in the following, the same symbol “cr” will be used to denote a cell assembly and the 
simplex that represents it) 071 {38]. The entire network can then be represented by a purely combinatorial 
simplicial complex T OH [32] whose maximal simplexes correspond to place cell assemblies 15911 , Sim- 
plexes in T may overlap: physiological studies demonstrate that a given place cell may be a part of many 
cell assemblies fiOl l4li . Many authors have suggested that place cell assemblies should overlap signifi¬ 
cantly in order to better represent contiguous spatial locations ll46l]49l : the more cells shared by cri and 0 - 2 , 
the closer the encoded locations are to one another. The most detailed representation of the environment is 
produced by a population of maximally overlapping cell assemblies, which differ by a single cell. In such 
case, a transition of the activity from one cell assembly cr\ to another crj occurs when one place cell in cr\ 
turns off and another cell in the new assembly o "2 turns on. The resulting simplicial complex T has the 
structure of a combinatorial d-dimensional simplicial manifold (in the literature also referred to as “pure 
complex” or “pseudomanifold”, OTTl ). 

Population activity in the cell assembly complex T. The simplicial complex T is a convenient instru¬ 
ment for relating place cell coactivity to the topology of the rat’s environment 1571 [38 11. The rat’s move¬ 
ments in the physical environment induce a packet of place cell activity that propagates in the hippocampal 
network—an “activity bump” E71 . In our model the propagation of the activity bump corresponds to an 
“active simplex” propagating through T. The resulting population activity vector is then 


~ C/cr, f<T,Vj 0 > fcr,Vi l > •••! fcr,Vi k )> 


CD 


where the first component f a represents the spiking rate of the readout neuron, and f T>V: denotes the spiking 
rate of the place cell c,- within the assembly cr. Roughly speaking, f^ Vi can be viewed as the firing rate of 
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Ci at the location where the place fields of the cells constituting the assembly cr overlap, which we refer 
to as the cell assembly field, l a (the domain /123 on Figure [l]3). A given place cell a is a part of many 
cell assemblies cri,cr 2 ,whose fields l (T] , l (Tl ,... are contained in the c,’s place field; thus, the higher the 
orders of the cell assemblies, the (statistically) smaller the /, r s Il39ll . Since the individual place cell spiking 
rates are well approximated by smooth Gaussian functions of the rat’s coordinates f36l . the quantities f a , Vi 
remain approximately constant over ! ir . The components of the population activity vector 0 in a given cell 
assembly can then be related to the corresponding place cells’ maximal firing rates f v . by a set of factors 
0 < h (TVi < 1 , that are specific to a given cell and a given cell assembly, 

/<r,v; = ha-, Vi f Vi , (2) 

which may be viewed as measures of the separation between the location l lT and the respective place fields’ 
centers (Figure |TJ3). In other words, the coefficients li lT ,v, provide a discrete description of the place field 
map’s geometry. 

As the rat moves from one cell assembly field to another (e.g., from l a - l to l (T ,,), the activity packet in T 
shifts from the maximal simplex cri to the maximal simplex cm, then to 0 - 3 , and so on, tracing a “simplicial 
path,” 


T = < 0 - 1 , cr 2 , (3) 

(these are “thick paths” in the terminology of f30ll ). As a result, every path y in the rat’s physical environ¬ 
ment corresponds to a simplicial path T 6 T, which can be viewed as an abstract representation of the place 
cell trajectory code used in EH . Flowever, in order to represent the path T in the hippocampal network, 
the activity of each place cell assembly cr e T should activate the corresponding readout neuron c (T . Thus, 
during the activation period, the net input from the presynaptic cells in cr should exceed the corresponding 
readout neuron’s firing threshold 6 a , 

^ w cr,Vif(T,Vi ^ 6a, (4) 

V/Scr 

where the coefficients w a , Vi represent the strengths of synaptic connections between the place cells and the 
readout neuron If 20 ll 6 lTl . 

In other words, this is a rate model in that the activity of cells is described by a single parameter: the 
tiring rate, /, related via coefficients li to the maximal rate ([2]). If the network is trained —the synaptic 
architecture is fixed, place fields are stable—then each cell assembly fires when the rat visits (or replays) a 
specific spot where the respective place fields overlap. Because this spot is very small compared to the size 
of place fields, the left side of 0 is the essentially the same every time. 

Using ([2]), the condition Q becomes 

^ b(r,Vifvi ^ @cr, (5) 

V/GfT 

where the variables 6 and 


b<r.v, — n'fr.v; /^rr.v, ? (6) 

are defined on all simplexes (i.e., all cell assemblies), and the variable / on the vertexes (i.e., place cell). 

Dressed cell assembly complex 7~*. The coefficients b fr , Vl can be regarded as characteristics of the 
maximal simplexes of T and the values f Vj as characteristics of its vertexes. Together, these parameters 
produce a “dressing” of the cell assembly complex with physiological information about the cells’ spiking 
and the network’s synaptic architecture. Equation 0 singles out a set $ 7 - of valid dressings, bj-, = {ba,v ■ 
eq. ([5]) is satisfied} which enable readout neurons to respond to presynaptic activity and thus defines the 
scope of working synaptic architectures of the place cell assembly networks. 
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FIG. 2 : Propagation of place cell activity along simplicial paths. A. The firing rate of the cell C2, required to 
activate the readout neuron of the cell assembly a~i, can be inferred from the firing rates of the other two cells, co and 
c i using the equation (JTOj) (red arrow over cr\). The the resulting rate of C2 and the rate /o of cq define the rate 
of the cell C3 required to activate the readout neuron in the next cell assembly 1x2 according to the equation (111 (red 
arrow over the oq). B. A network of maximally overlapping cell assemblies is represented by a simplicial manifold 
T t . The replayed sequences correspond to simplicial paths that can be closed (T1) or open ^2). 


Replays occur on a millisecond time scale and produce only a few spikes per activity period llT9il (Fig¬ 
ure Gp). which is comparable to the stochastic background activity of neurons. In order to distinguish 
cell assembly activation from the assembly’s background activity, the readout neuron should be tuned to a 
particular combination of inputs; the physiological process most likely involves gating specific parts of the 
dendritic tree by timed inputs from the presynaptic cells. Here, the model employs a simplified version of 
this process: we propose that the readout neuron should remain in a sensitive, near-threshold state ll50l that 
allows it to quickly respond to the cell assembly during each individual step of replay. Thus, we hypothesize 
that during replays the inequality ([5]) may be approximated by the equation 

2 V,v,/v, = e<r, (7) 

V/GCT 

which further restricts the set of valid dressings to a special marginal set S 7 -,, for which spontaneous replays 
of place cell assemblies are possible. Note, however, that equation 0 does not fix the values of the param¬ 
eters f v . and bo-^, and so it allows a significant variability of spiking activity and of the synaptic connection 
strengths. To emphasize the assumption that place cell activity during replays elicits a response from the 
readout neuron at a constant rate fo- > 0 , it is convenient to use the notation 

d,r = ba-'O-fa, ( 8 ) 

where b^a- > 0 is a fixed parameter that can be interpreted as the readout neuron’s susceptibility to dis¬ 
charge. 

Replays. Equation ([7]) defined at each simplex of T, lf30l provides a simple tool for building a model 
of hippocampal replay. As an illustration, consider the case when 7~* is two-dimensional and let cr\ = 
[vo, vi, V 2 ] be a 2D simplex representing an assembly of three cells with the firing rates f, 0 , f Vl and f V2 . If 
equation ([7]) holds over cr \, then the readout neuron c (Tx fires with the rate f T] in response to the coactivity 
of ci, C 2 and C 3 , 


b(Tl,Vofvo + bo- 1; V1 


f 


vi b (T ^ V2 fv2 — bo-i,crifiTi- 


(9) 


Suppose that equation (|7} also holds for an adjacent (maximally overlapping) cell assembly, represented by 
an adjacent simplex 0-2 = [vq, V 2 , V 3 ], so that the second readout neuron c (r , fires with the rate f Tl 


bfr-.V] f\ : ] T b a - 2 ,Vlf\ , 2 b(T2,V 2 fv3 ~ b(j-2,(T 2 f(T2 • 


( 10 ) 
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A key observation here is that, since <X 2 shares vertexes v\ and V 2 with cri, the corresponding firing rates 
f v j and f V2 in ( p~Q] > define uniquely the firing rate of the remaining cell, / V3 , required to activate the readout 
neuron c a2 (Fig.[2}\), 


fv 3 


y O- 2 ,V 3 


~(b(T 2,C2 fo~2 b(T2 ,V 1 /v 1 ' l ? cr 2 , V2 f> 2 )• 


( 11 ) 


Similarly, if there is another 2D simplex cr 3 = [V 2 ,V 3 ,V 4 ] adjacent to o" 2 , then, once the value /,, 3 is found 
from ( [IT] ), the firing rate at V 4 can be obtained from f V2 and f V} , and so on (Figure [2j3). In other words, 
once the synaptic connections b& tV are specified for all simplexes, equation ([7]) can be used to describe the 
conditions for transferring the activity vector f (r over the entire complex 7~* |30| . Notice however, that 
equations (|9j)-(pT|) do not specify the mechanism responsible for generating place cell activity; they only 
describe the conditions required to ignite the cell assemblies in a particular sequence. While the subsequent 
simplexes cr n and cr n+ 1 in the simplicial path ([3]) are not necessarily adjacent, the activity according to 
equation (|7]> is propagated along a sequence of adjacent maximal simplexes, such as depicted in Figure [2^3. 

Discrete Holonomy. Using the notation 




CT 

x,y 


b<j,x 

b<r,y 


equation ([IT]) defined over a simplex cr p can be rewritten in matrix form 


f = M r,Vs f 

L q — 1Y1 q,p 1 / ? ’ 


( 12 ) 


(13) 


where the “transfer matrix” M q v p propagates the population activity vector from the incoming facet of the 
simplex cr p into the activity vector of the outgoing, opposite facet shared with the next simplex <j q (edges 
[vo,vi] and [vj,V 2 ] respectively on Figure [2j\), in which the vertex y s of the simplex cr p shuts off and the 
vertex v t of the adjacent simplex activates. If there is a total of n simplexes in the path F (n = 17 for 
the closed simplicial path T\ shown on Figure |2j3) then the corresponding chain of n equations ( [13] ) will 
produce 




(14) 


If the simplicial path is closed, then the activity vector should be restored upon completing the loop, i.e., 
f n = fi. According to ( |T4] ). this will happen if the product of the transfer matrices along T yields a unit 
matrix, 


= M Vl " Vs " M V, "~ 


(15) 


It can be directly verified, however, that condition ( fl5| ) is not satisfied automatically: the product of transfer 
matrices (T3j) has the structure 


M r = 


10 O' 
0 1 0 

/era kt, 2 1 + Kr,3 , 


(16) 


which differs from the unit matrix M\ # 1 (see Appendix). This implies that a population activity vector 
icr is in general altered by translations around closed simplicial paths, f,Aktau) + U-(ktnd)- To formulate 
this another way, the spiking condition ([7]) does not automatically guarantee that the readout neurons will 
consistently represent spatial connectivity; the latter requires additional constraints (151, irrespective of the 
mechanism that shifts the activity bump. 

Mathematically, a mismatch between the starting and the ending orientation of the population activity 
vector is akin to the differential-geometric notion of holonomy which, on Riemannian manifolds, measures 
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B. C. 



FIG. 3: Discrete geometry of a dressed simplicial complex. A. Discrete holonomy: a population activity vector 
(red arrow) changes its direction from simplex to simplex as described by ( fl3j ). Upon completing a closed path, the 
starting and ending vectors may differ, f slar , + which indicates nonzero holonomy. B. A 2D elementary closed 
path of the order n encircling a vertex vq. The “pivot” vertex Vo carries the discrete curvature coefficients defined by 
0- c. A higher dimensional elementary closed path consisting of ^/-dimensional simplexes (one such exemplary 
simplex cr 1 ^ is shadowed) sharing the same (d - 2)-dimensional face, the pivot simplex <ri, shown in red. The d — 2 
dimensional pivot simplex <x shown in red carries the curvature coefficients K,- r i . 


the change of a vector’s orientation as a result of a parallel transport around a closed loop If30l I33l435ll . 
Hence, the requirement ( p~5| ) that the activity vector should be the same after completing a closed simplicial 
trajectory implies that the discrete holonomy along paths in 7~* should vanish. 

Discrete Curvature. In differential geometry, zero holonomy on a Riemannian manifold is achieved 
by requiring that the Riemannian curvature tensor R‘. k/ associated with the connection P.^ vanishes at every 
point v EH351. This condition is established by contracting closed paths to infinitesimally small loops 
encircling a point x and translating in parallel a unit vector ft around that loop. The difference between the 
starting and the ending orientations of ft defines the curvature at the point x If35l . An analogous procedure 
can be performed on a discrete manifold 7~*. However, there is a natural limit to shrinking simplicial paths: 
in a (/-dimensional complex, the tightest simplicial paths consist of (/-dimensional simplexes cr (d) which 
intersect the same (d - 2) dimensional face (see Figure |3j3). Such a path T^d- 2 ) we will call an “elementary 
closed path”, following l l30l l. The order s n of such a path is defined by the number of (/-dimensional 
simplexes cr {d) encircling a simplex cr w ~ 2) . In the following we will use the short notation & for the “pivot” 
simplexes & = cr i ' d ~ T> whereas the elementary simplicial path encircling & will be denoted as p r . 

In order to ensure zero holonomy of place cell activity along all closed paths in 7~*, it is sufficient to 
verify that the holonomy vanishes for all elementary closed paths l(30ll . The product of the matrices 
encircling the pivot <x (Figure |3j3) has the same form as equation ( [T6| ); however, the coefficients k,- t j at the 
bottom row of the matrix M,- r can be viewed as the curvatures defined at <x. Thus, to ensure zero holonomies, 
the conditions 


kx,i — 0, 


(17) 


i — 1, ..., <7+1, must be imposed on the connection coefficients b, jyi at every pivot simplex <x of a <7- 
dimensional dressed cell assembly simplicial complex. For example, an elementary 2D closed path encir¬ 
cling a vertex vq with n simplexes enumerated as shown on Figure [3]7 yields the holonomy matrix 


AT, 


vo 


'10 O' 
0 1 0 

, A’o.l G’<>,2 1 + Kvo,3 „ 


The values k V() j, i = 1, 2, 3, of the bottom row that distinguish M Vo from the unit matrix should be considered 
as discrete curvatures defined at the pivot vertex vo (see Figure [3p and l l30l ), which need to vanish in order 
to ensure a consistent representation of space during replays. 





















Since there exists a finite number of pivot simplexes, the number of constraints <0 on a given dressing 
bj- e St. is finite. Thus, the scope of nontrivial zero holonomy conditions ( p~5| ) drastically reduces and 
the task of ensuring consistency of translations of the population activity vectors over T* becomes tractable. 
Nevertheless, zero curvature conditions ( fT7j ) are in general quite restrictive and impose nontrivial constraints 
on the synaptic architecture of the place cell assemblies. As the simplest illustration, consider the case 
when the firing rates of all the place cells and readout neurons are the same: f Vi = f (Jk = /, and all 
the connection strengths from the place cells to the readout neuron in all cell assemblies are identical: 
= b, giving a constant connection dressing Bt, ■ It can be shown that in this case the resulting transfer 


'm.V; 


matrix is idempotent, that is M~ = 1, so that the zero curvature condition (171 is satisfied identically for the 
even order elementary closed paths and cannot be satisfied if the paths’ order is odd. Under more general 
and physiologically more plausible assumptions equation ( |T7| ) does not necessarily restrict the order of the 
cell assemblies. However, the domain of permissible dressings, St . is significantly restricted by (17 1 , as 
compared to the domain occupied by the synaptic parameters of the unconstrained cell assembly networks. 


III. STATISTICS OF SYNAPTIC WEIGHTS IN THE LIMIT OF WEAK SYNAPTIC NOISE 


The zero curvature constraints ( [17] ) affect the net statistics of the synaptic weights. Since the structure 
of the full space of marginal dressings Sr. and of the corresponding probability measures is too complex, 
we considered a family of connections parametrized as 

b(T,v = b (l + E a , v + , (IS) 


in which the fluctuations s ay are normally distributed, 

1 4-,v 

P(£<r,v) = — -p e fl ■ 

£ \7T 


(19) 


In the absence of zero curvature constraints, cell assemblies are uncoupled and the synaptic fluctuations are 
statistically independent, so that the joint probability distribution of s ay is 


ne) = Y\ P ^- ( 2 °) 

<T,V 

Under zero-curvature conditions ( [T7| ) the parameters of the synaptic architecture are coupled (Figure [5]) and 
the probability distribution for a particular variable £, ry is obtained by averaging the joint distribution ( |20| ) 
under delta-constraints: 


/y-% d~\~ 1 

■J \^\^6(k & /£))P(£)D'£, 


& i= 1 


( 21 ) 


where C is the normalization constant and D's = ds^ denotes integration over all £ <T >y + £ a ,v 

In the Appendix it is demonstrated that for weak fluctuations, the shape of the distribution ( |T9j ) remains 
Gaussian, 


P ( e cr,v) — 



( 22 ) 


but its width decreases: £* < e. Thus, zero curvature conditions narrow the distribution of the uncorrelated 
weights, i.e., produce a “tuning” of the synaptic connections b tjy . This result also applies to the synaptic 
weights: in cases where the place fields are distributed regularly, so that the coefficients h ljyi have a well 
defined mean, h , and a small multiplicative variance, 


b(T. V, — h (1 + b,T.\ ) , 
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5cr,Vj c 1, the coefficients fj^ V2 are approximately defined by the ratios of the synaptic weights, 




CT 

Vl,V 2 


b(T,V 1 


Wcr,vi htr.V] 
W(T,V2 h<T,V2 


W (T,V l /, c r \ 

(1 + 0 (j- jVl — 

n cr,V2 


Wcr.vi 

^<T,V 2 


and therefore the zero curvature conditions produce the same effect on w (T ^, as on b (J . v , i.e., reduce the 
variability of synaptic weights. 

Understanding the effects produced by zero curvature constraints ( fF7| ) on a wider range of fluctuations 
is mathematically more challenging. The qualitative results obtained here, however, may generalize be¬ 
yond the limit of small multiplicative synaptic noise and could eventually be experimentally verified. A 
physiological implication of the result (22) is that the distribution of the unconstrained synaptic weights in 
a network that does not encode a representation of space (e.g., measured in vitro ) should be broader than 
the distribution measured in vivo in healthy animals, which can be tested once such measurements become 
technically possible. 


IV. DISCUSSION 

The task of encoding a consistent map of the environment imposes a system of constraints on the hip¬ 
pocampal network (i.e., on the coefficients b lT ^) that enforce the correspondence between place cell activity 
and the animal’s location in the physical world. Here we show that zero holonomy is a key condition, which 
is implemented by requiring that curvatures vanish at the pivot simplexes. This approach works within a 
combinatorial framework, but a similar intuition guided a geometric approach 1291 . where the place cells’ 
ability to encode the location of the animal—but not the path leading to that location—was achieved by im¬ 
posing the conditions of Stake’s theorem If34l on the synaptic weights of the hippocampal network, which 
were viewed as functions of Cartesian coordinates. Our model is based on the same requirement of path- 
invariance of place cell population activity, implemented on a discrete representation of space—a dressed 
abstract simplicial complex T t —without involving geometric information about the animal’s environment. 

In particular, note that the concepts of “curvature” and “holonomy” are defined in combinatorial, not 
geometric, terms. This is an advantage in light of (and indeed was motivated by) recent work indicating that 
the hippocampus provides a topological framework for spatial experience rather than Cartesian map of the 
environment B4L and it also makes our model somewhat more realistic. It does, however, lead to a number 
of technical complications. For example, discrete connections ([6]) defined over 7~* are nonabelian PUL so 
using the approach of ||29l would require a nontrivial generalization of Stake’s theorem, which is valid only 
in spaces with abelian differential-geometric connections l5ll . Our approach is based on the analysis of 
discrete holonomies suggested in the pioneering work of POl which, in fact, explains the mathematical 
underpinning of the Stoke’s theorem approach in both abelian and nonabelian cases P3I4351 . Indeed, the 
zero-holonomy constraint ensures that no matter what direction the activity is propagated in the network 
(forward, backward, or skipping over some cell assemblies), the integrity of the spatial information remains 
intact. 

Generality of the approach. A key instrument of our analyses is equation ([7]), which describes the 
conditions necessary for propagating spiking conditions over the cell assembly network. The exact form 
of this equation is not essential—a physiologically more detailed description of near-threshold neuronal 
spiking lf52f - l54ll could be used to establish more accurate zero holonomy and curvature constraints on the 
hippocampal network’s synaptic architecture, which should be viewed as a general requirement for any 
spatial replay model. 

The assumption of maximally-overlapping place cell assemblies may also be relaxed, since equation (|7|) 
can be applied in cases where the order of the cell assemblies varies, that is, when the simplicial complex 
7~* is not a manifold but a quasimanifold (see Figure [4] and li55l l56l ). Unfortunately, implementing the 
“zero holonomy” principle in this case would require rather arduous combinatorial analysis. For example, 
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FIG. 4: A replay in simplicial quasimanifold. An example of a simplicial quasi-manifold containing 2D and 3D 
simplexes. The activity of cells in the 3D simplexes is induced from the 2D simplexes approaching its sides. Two 
simplicial paths are shown by gray triangles, marked by red dotted lines. 


propagating the activity packets using ( [Ft] ) would impose relationships between the dimensionalities of the 
maximal simplexes and their placement in 7"*, i.e., require a particular cell assembly network architecture. 

Learning the constraints. In this paper, the requirements ( |T7] ) and (17 1 enforcing path consistency of 
place cell replay are imposed on a fully trained network: it is assumed that the place fields have had time 
to stabilize |[59l and that the cell assemblies with constant weights w a - !Vj have had time to form |j4[. In a 
more realistic approach, these constraints should modulate the hippocampal network’s training process. For 
example, if the unconstrained network is trained by minimizing a certain cost functional S(bj -,) Ifl8l 1621 
then the constraints ([17]) would contribute an additional “curvature term” Ribq-J, 


S(b T ,)-> S(b r ,) + R(b Tt ), 
defined, e.g., via Lagrange multipliers r,-, 

R(brJ = = Y 


r; Ka 


(23) 


(24) 


&eT, 


Physiologically, the network may be trained by “ringing out” the violations of the conditions ( [17] ) in 
the neuronal circuit, i.e., by replaying sequences and adapting the synaptic weights to get rid of the centers 


of non-vanishing holonomy. Curiously, the role played by R(bj-J in (23) resembles the role played by the 


curvature term in the Flilbert-Einstein action of General Relativity Theory | |34l . which ensures that, in the 
absence of gravitational field sources, the solution of the Flilbert-Einstein equations describes a flat space- 
time. By analogy, the constraints imposed by ( [17] ) may be viewed as conditions that enforce “synaptic 
flatness” of the hippocampal cognitive map. 

It is worth noting that the mechanism suggested here is an implementation of the zero holonomy condi¬ 
tion in this simplest case of the reader-centric cell assembly theory that is consistent with physiology. The 
place cell readout might involve, instead of a single neuron, a small network of a few neurons (not yet identi¬ 
fied experimentally), which might require a different implementation of zero holonomy principle, depending 
on the specific architecture of such a network. If the readout network is a cluster of synchronously activated 
downstream neurons, then this cluster of cells could be viewed as a “meta-neuron” and the proposed ap¬ 
proach would apply to this case as well. More complicated architectures would require modifications, but 
it is reasonable that the reproducibility of the population vector would require zero holonomy in all cases. 
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VI. APPENDIX 


Transfer matrix construction is carried out for the 2D case, since higher dimensions arc similar. In the 


matrix form, equation (111 defined over the simplex cr p , can be written as i q = M qp J f p , in which the matrix 


M VkV] = 

lvl q,p 


0 

1 


o 
0 


Fo>,n- ^vi,vk A* VuVk ) 


Vj,Vk 


transfers the activity vector 


f T 

p 


(fo-p, fvi, fvj) , 


defined over the incoming edge [v,-, vj\ of the p-th simplex into the activity vector at the outgoing edge 
[v;, V*] of the same simplex (e.g., from the edge [vq, vi ] to the edge [vi, vi\ of cr\ on Figure |2]A), 


f T 

p 


(UpJ Vi J Vk )- (25) 

To ignite the readout neuron of the next cell assembly cr q , which shares the edge [v ? , vy. ] with cr p the vector 


(25 i needs to be transformed into 


f T 

9 


fvil f\’ k ) 


by the diagonal matrix D pq = diag(/) r / f a ,1,1). Together, these two operations produce the transfer 
matrix 


M nVi = D M nVj = 
lyl q,p LJ P,q lvl q,p 


fo- q lfo- p 0 


0 ' 

0 1 0 

r p _ O-p _ Up 

H(T p ,V k /U.-.V* f*Vj,V k 


A direct verification shows that a product of n transfer matrices that stat and end at the same simplex, has 
the form ( p~6] >, in which /q-, are nth order polynomials of the coefficients (12 1 . 

Tuning of the fluctuation distribution. In case when the fluctuations are small, e <s: 1, the constraints 
(fTT]) uncouple (Figure [5]) yielding linearized curvature coefficients 


Kfr,i — 


Q&.V,(T, 


y Scr,v, 


(26) 


vstA 


where f2oy ';< T , v arc constant coefficients and the summation is over the vertexes of the elementary path V,-, 
(Figure [3j3). To simplify the expression ([20]), we rewrite it using indexes p = (cr, v) and / = (<x, i), 


Hsp) = C f • ■ ■ f PI SUiiSpPyme,/) n ds p ', 

J J i p'*p 


and exponentiate the delta-functions, 


OO 

P{s p ) = C f-f f -f ns,/) P| drpe™ P| ds,y. 

J J q=-J ’ p'*p 


( 27 ) 


Using the joint distribution (20 1 and the linearized expressions (26 1 in ([27]) produces 


P(s p ) = Ce 


OO 

J-/n drne^p /.../n* 

T]=-O0 l TT (U _ \ P 


ds P ’e ^ e l1 > lV >p’ s p\ 


u e (br,) 
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FIG. 5: Coupling between simplexes. A schematic illustration of a maximal simplex cr\ span by three pivot vertexes, 
vo, v\ and V 2 , and shared by three overlapping elementary paths F [, To and F 3 in a 2D cell assembly complex. As a 
result, the synaptic connectivity coefficients ba- liV will appear in three sets of discrete curvatures, k VI]A , k Vi j and k V 2 j, 
i = 1, 2, 3, bootstrapping the constraints (| 1 1\. 


where the Vj p are the coefficients obtained by collecting the terms proportional to e p >s produced by (26). 
Completing the square and integrating over s P ' yields a Gaussian integral over a positive quadratic form 
A = EE T , 


P(s p ) = e 


, °° 

f ,'l I 

rj=-co 1 


drne^We-** 2 ^. 


where v p is the pth row of the matrix V. Evaluating (28) yields 


P(s p ) = Ce E 


JL _ LL V A~ l v 


= Ce 


(28) 


(29) 


where 


2 2 
£, = £ 


(1 


+ v p A~ 




(30) 


Since the second term in the parentheses is positive, £* < £, which indicates narrowing of the uncoupled 
distribution ( fT9| ). The magnitude of the correction in ( |30| ) depends on the topological structure of the 
coactivity complex (e.g., its dimensionality d and the statistics of the pivots’ orders, n) and on the dressing 
parameters, 'Bq-. In the approximation (18), £ <sc 1, the diaginal matrix elemens of the matrix A are of the 
order dh, and hence the v p A~ l v p ~ dIn. 
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